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In the simplistic quantum mechanical picture of flavor mixing, conditions on the maximum size 
and minimum coherence time of the source and detector regions for the observation of interference — 
as well as the very viability of the approach — can only be argued in an ad hoc way from principles 
external to the formalism itself. To examine these conditions in a more fundamental way, the 
quantum field theoretical S-matrix approach is employed in this paper, without the unreahstic 
assumption of microscopic stationarity. The fully normalized, time-dependent neutrino flavor mixing 
event rates presented here automatically reveal the coherence conditions in a natural, self-contained, 
and physically unambiguous way, while quantitatively describing the transition to their failure. 
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Q ■ I. INTRODUCTION 

m : 

psj ' Several recent works have examined neutrino flavor mixing by considering the neutrino production/mixing/ detection 
T— I , as a single process in the context of quantum field theory (QFT) ||l|-^ . Such a framework clarifies several conceptual 
difficulties associated with the familiar quantum mechanical (QM) model of the flavor mixing process (see Ref. |£| for 
a listing of some of these) . One conceptual difficulty associated with the simphfied QM picture is that it postulates 

• neutrino flavor eigenstates of indefinite mass, while in QFT external particles (asymptotic states) are generally required 
\ to be on-shell. Hence the usual methods of calculating neutrino production rates in QFT would be rates for neutrinos 

^f-^ • of a particular mass, precluding interference between neutrino states of different mass. In a QFT description of a 
^\ ' neutrino mixing experiment, this problem is resolved by considering the neutrinos to be virtual particles. After all, it is 
I the measurable, on-shell external charged leptons associated with the neutrino production and detection processes that 
0^ operationally define what is meant by "neutrino flavor mixing" ; the neutrinos themselves are not directly observed. In 
the relativistic limit, the same factors that constitute the "oscillation amplitude" in the simplified quantum mechanical 
picture can be identified in the amplitude for the overall neutrino production/mixing/detection process. 

While descriptions of neutrino fiavor mixing in QFT have provided insight, some shortcomings remain. As noted 
in Ref. one problem is that the calculations ffj-Ql are not carried out to normalized event rates. Without normal- 
^ ization, one cannot definitely say that one has identified an "oscillation probability. "|^ In addition, Refs. are 
restricted to particular neutrino production and detection reactions, while Ref. [Q employs idealized two-state systems 
^ ■ as source and detector. Since one would hope to justify the use of the simple QM model in general circumstances, 
. , such restrictions should not be required. 

• Another potential pitfall is a failure to distinguish between macroscopic stationarity and microscopic stationarity.^ 
Some previous studies invoke microscopic stationarity, either implicitly or with explicit reference to bound 
states 1^,^ in the source and/or detector. While sources and detectors as a whole can sometimes be considered 
stationary on a macroscopic basis, the claim that individualparticles in the source and/or detector remain unperturbed 
in coherent states over macroscopic time scales is dubious.^ A good example is the Sun: While macroscopic variables 
such as density, pressure, and so on may be stationary, zooming in to atomic scales one sees a roiling thermodynamic 
bath of particles being created, destroyed, and scattered on rapid time scales. Clearly there is no hope of appealing 
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^A normalized probability is given in Ref. but it is not the experimentally relevant one, which is identified by making a com- 
plete connection of the squared amplitude for the production/mixing/detection process with the form (neutrino flux) (oscillation 
probability) (neutrino cross section) [^. 

^This critique applies only to QFT analyses of mixing experiments, which by their very nature purport to describe microscopic 
processes in the source and detector. Once the feasibility of the use of flavor eigenstates in a simplifled model is established, 
macroscopic stationarity can be sensibly employed, as in Ref. | ]lo{ . 

^The "macroscopic time scale" at issue here is the signal travel time between source and detector. In the case of astrophysical 
neutrino sources, this time is macroscopic indeed. 
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to bound states of the nuclei and electrons which collide to produce neutrinos in the Sun! Even in detectors which 
have bound state particles, it is difficult to conceive of these states as being coherent over macroscopic time scales. 
For example, water Cerenkov detectors see charged lepton wave packets with finite energy and time spread (after all, 
an "event time" is recorded whose uncertainty is much smaller than, for example, the signal travel time between the 
Sun and Earth). While the overall chain of the processes of detection can be rather complex, one would expect that 
at least part of the reason for the finite spread of these detected charged leptons is the limited coherence time of the 
bound state particles with which the initial detection interaction takes place. 

In short, if one employs a QFT description of flavor oscillations in order to overcome the conceptual difficulties of 
the simplified QM model, one should also pay the price of being realistic about the lack of microscopic stationarity 
in order to complete a convincing picture. In this paper, this is accomplished by treating all of the initial and final 
state external particles as wave packets which have finite overlap in space and time in the source and detector. Thus 
it is similar in spirit to Refs. but in pursuit of generality the treatment does not specify particular neutrino 

production and detection mechanisms or specific functional forms for the wave packets of the external particles involved 
in the production and detection processes. Another difference of the present treatment is a greater emphasis on the 
coordinate space Green's function, as in Ref. where it was used to make direct contact with the standard simple 
coordinate space formalism for the MSW effect. In this work no integration is performed over the time coordinate 
of the detection event, since many experiments — including those employing water Cerenkov detectors — record this 
time. (An integration is performed over the unobserved time coordinate of the emission event.) Finally, the detailed 
connection of the squared amplitude for the microscopic neutrino production/mixing/detection process to macroscopic 
event rates will be made. Only this complete connection — with all factors accounted for — enables one to define an 
oscillation probability. 

II. S-MATRIX APPROACH TO NEUTRINO MIXING PROCESSES 

A common application of QFT is the description of particle collisions in accelerators. Rates or cross sections 
associated with these processes can be obtained in a heuristic manner directly from the plane wave scattering S- 
matrix computed from Feynman diagrams. 



S{{p}) - 1 EE (2^)4^4 Y^i-l f^p, ^Map}). 



(1) 



In this expression {p} is the set of external particle momenta pi, di = \ for incoming and for outgoing particles, and 
M. is the 6 function- free matrix element. The event rate obtained from Eq. (nl) isF 



dV = {2^f5' V'-'\M{{p})\ 



n 



1 



2£;n 



n 



dp,; 



(27r)3(2£;p^,) 



(2) 



Here V is the three- volume in which the entire process occurs, / and F are the numbers of particles in the initial and 
final states, and the components of the on-shell four- moment a are (iJp. ,pj). This mnemonic for arriving at event 
rates is possible because the interactions of interest occur in a single, small spacetime volume. It is more convincingly 
justified, however, by a wave packet description (e.g., Ref. fl^). 

One of the reasons one considers a QFT description of neutrino flavor mixing is that the standard picture of 
requiring external particles (asymptotic states) to be on-shell precludes the existence of massive neutrino flavor 
eigenstates (assuming that each charged lepton couples to multiple neutrino fields of different masses). Accordingly, 
one considers the neutrinos as virtual particles in a Feynman diagram in which the charged leptons at the source 
and detection vertices identify the neutrino flavor. Neutrino flavor mixing then results from interference of diagrams 
whose intermediate neutrinos have different masses. 

In this picture it is not possible to compute event rates directly from the S-matrix with the usual mnemonic 
described above. This is because a neutrino oscillation experiment involves neutrino production and detection regions 
which are widely separated in space. In contrast to the case of accelerator particle collisions, the interactions of 
interest do not all occur in a single volume element. In addition, as argued in Sec. ||, in this microscopic picture the 
production and detection of a single neutrino will be separated in time as well as space. 



*The conventions for the metric, 7 matrices, and normalizations employed here are the same as those of Ref. [ pA[ . 
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In order to describe the spacetime localization one must fall back on a wave packet description of the external 
particles, in which the amplitude is a superposition of plane wave amplitudes: 

„1+Fd Is+Fs 

^ = jll idPj] (bojiPj^P^) n [dH^s^ikM [^({fc},M)-l], (3) 

j i 

where (for example) [dpj] = rfpj/ [(27r)'^-^2ii^pJ , {k} are the external momenta connected to the vertex causing 
neutrino production, and {p\ are the external momenta connected to the vertex associated with neutrino detection. 
The quantities {k} and {p} denote the peak of the wave packets' distribution of three-momenta. There are Is 
incoming and Fs outgoing particles connected to the production vertex, and 1 incoming and Fd outgoing external 
particles at the detection vertex. The origin of the source wave packets is taken to be the spacetime point xs'- 

(j^S^^K^d = as,(k, -k,)e-*(-l)'"^--^^ (4) 

and similarly for the detector wave packets with origin yo- The real function a(k — k) is peaked about k. In order 
that Eq. describe the amplitude for interaction of one localized external particle of each type, the wave packet 
normalization must be, e.g. for the source packets, 

-^^\4>si{h,'ki)f' = j -^^\asi{^i-ki)? = 1- (5) 
It is convenient to define the transform 

V'k.(x) = j ^a5.(k, -kJe-(-i)"''-, (6) 

where Si behaves the same as the di of Eq. (^. The normalization of this function is 

ydx|V'k,(x)i'=i, (7) 

which follows from Eq. (||). 

The next step is to transform the momentum-based Eq. into a coordinate space expression. The S'-matrix in 
Eq. (0) can be expressed 



S{{k}, M) -1 = J d^y e'E,(-i)''P' S' J d'^x E,(-i)°"='- 



s 



J ■^^^'"'^'''"^^'I^PlG{s)PrMi, (8) 

in which s is the off-shell neutrino propagator momentum. The upper (lower) sign of =F in the exponential is 
for neutrino (antineutrino) mixing. This arises from choosing x (y) to always correspond to the source (detec- 
tor). That is, for neutrino oscillations of flavor a to flavor /3, the propagator is iG^"{y,x) — {T{i'^{y)D"{x)})Q = 
i J d'^s {2'k)~'^ e~^^'^y~^'>G'^°'{s) (with T{} and ()o denoting a time-ordered product and vacuum expectation value 
respectively), while for antineutrino oscillations a (3, the labeling is iG"^ {x,y). V — A interactions have been 
assumed; Pl and P^i are the left- and right-handed projection operators, with Mi and M2 column and row vectors in 
spinor space respectively. For neutrino mixing, Mi — Afi({fc}) and M2 — M2{{p}) are respectively associated with 
the neutrino production and detection reactions. For antineutrino mixing, M2 = M2{{k}) and Mi = Mi{{p}) are 
respectively associated with the production and detection reactions. The partially transformed propagator G'(s°, y, x) 
is defined by 

■ e^''<y~^^G{s) = / ^ e^''°(y°-^"'> G(s°, y, x). (9) 



(27r)4 ^ ' J 2t: 

It is assumed that the wave packets a(k — k) arc sufficiently well-peaked that integrals of the following form can be 
evaluated in an approximate manner: 



/ 



[rffc.]</.5.(fc„k,)e'(-i)"'='-M(k,) = / as.(k, -kJe*(-i)"'=-(----)M(kO 

V-k, ((x - X5) - {x° - 4)vk j Af(k,). (10) 
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InEq. @, Vk. is the wave packet's group velocity (Vkifc°)|k,=k. = k,:/-E^k , and wave packet spreading has been 
neglected. Similar expressions hold for the detector wave packets. With this approximation, the amplitude of Eq. (||) 
becomes 



^Is+Fs 



1 + Fd 



2E„ 



'Is+Fs 

W V^k^ ((x - X5) - - )vk 



n ((y - yd) - [if - yl)^^^ 
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l + Fo 



(11) 



in which the bars in Mj^ and Af o signify that these quantities have been evaluated at the peak momenta of the wave 
packets. 



III. APPROXIMATION OF THE WAVE PACKET OVERLAP 



It is convenient at this stage to adopt an approximation regarding the overlap of the wave packets that captures 
the essential physics while maintaining mathematical simplicity. The initial and final state wave packets in the source 
(for example), traveling with their various group velocities, overlap in a limited region of space for a limited time. To 
give a specific definition to this spacetime volume of the overlap, Vs, centered on xg, it is convenient to define 



Esix - xs) = 



Fs+Fs 



n ^k.(x5,4) 



Fs+Fs 



n ^k. ((x - X5) - - 4)vkj 



where the notation 



V-k^ (xs, x%) = Vk, ((x - Xs) - - x%)-^^^ 
has been adopted in the first factor. Then Vg is defined by 



(12) 



(13) 



d'^x {Es(x ~ xs)\ = 



d X cxp <^ 2 In 
^Det [{Ws)^,.] 



l-^{Ws)^Ax~xsrix~xsr 



where 



dx>^dx'^ 



Esix - Xs) 



(14) 



(15) 



The timelike (spacelike) components of iyVs)^v refiect the spread of energy (momentum) available in the reaction, 
while the timelike (spacelike) components of iyVg^^y^i, characterize the extent in time (space) of the wave packet 
overlap. Similar considerations apply to the detector region. 

It is only necessary here to consider the Green's function for neutrino propagation through the vacuum. Inspection 
of Ref. 1^ indicates that the vacuum propagator results will be applicable in a relatively direct way to the case of 
neutrino propagation through a medium of constant density. Generalization to the case of a medium of varying density 
would be more complicated, however. While interference terms for this case have been calculated in the context of 
the simple quantum mechanical model [p^ , they are not relevant to current observations of astrophysical neutrinos. 
Hence the effort to study the microscopic origin for the damping of interference terms already deemed irrelevant does 
not seem to be worthwhile at present. 

Focusing on the vacuum case — for which the interference terms are of current experimental interest — the final 
factors in Eq. dll]) can be expressed as 



(16) 
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where and Mj^ are respectively the two-component subspinors that remain nonzero in Mj^Pl and PrM_^, and 
G(s°, y, x) is the nonzero 2x2 submatrix in PlG{s^, y, x)Pr. Because the overlaps of the wave packets are restricted 
to the vicinity of xg and yo, and because \yD — xs| Ls,Ld, the leading contribution from the Green's function 
is of the form [01 



G"'3(s0, y, x) « - ^ U^kU;, ( s° - ^ <T . L 



giSfcL-(y-x) 

4vr|yn - X5I 



(17) 



In this expression, the flavor and mass fields are related by Va{x) — '^f^Uak4'k{x); Sk — \/ (s^y^ — rn^, in which 
TOfc is the mass associated with the neutrino field ^k{x); cr is the three- vector of Pauli matrices; and the vector 
L = (yo — X5)/|y£) — X5I points from the source to the detector. For neutrino oscillations, /3 is the flavor of the 
charged lepton associated with the source reaction, and a is the flavor of the charged lepton associated with the 
detection reaction. For antineutrino oscillations these assignments are reversed. 

With these preparations the remaining integrations in Eq. ( [ll| ) can be performed. Employing Eqs. (12), ([l^), and 
(p^), and employing a similar approximation to that employed in Eq. ( |l^ ) for the x and y integrals, the amplitude 
becomes 



-'^^^ (X5,4) 



n 



2£;. 



(4Vs)(4Vg) 
4vr|y£, - xsl 



(18) 



in which the function cxp[— Z?fc(s'^)], with 



(19) 



enforces energy-momentum conservation to the extent allowed by the finite overlap in space and time of the external 
particle wave packets. In Eq. (p^), the notation 



PD 



I 

^+J2i-^f'Pr 



Cfc EE ±5°, L./(sO)2 



(20) 
(21) 

(22) 

has been employed. 

If plane wave final states and stationary initial source and/or detector particle states had been employed, the finite 
energy spread indicated in Eq. ( p^ would have been replaced by an energy delta function, suggesting the idea that 
the neutrinos are energy eigenstates. In such a case, one finds s° = ^P% = ±^5 (it will be recalled that the upper sign 
is for neutrino emission at the source, and the lower sign for antineutrino emission). It is easy to see by considering 
sample neutrino production and detection processes that kg and p'j^ should be positive quantities. 



While Eq. (19) indicates that a range of s contributes to the amplitude, there is a value of s — call it (s )/c — which 



makes the largest contribution to the amplitude. This is the value of for which Dk{s'^) has its minimum value. 
The relative degrees to which overall energy and momentum are conserved compete in determining (s°)fe. Since 
the external particles travel at speeds less than the speed of light, however, the timelike components of the tensors 
{Wg]-i)fj,i, will be larger than the spacelike components. This means that, in analogy with the stationary situation 
mentioned above, {s^)k may be taken to be positive (negative) for neutrino (antineutrino) emission at the source. 
This, together with the fact that only the region in the vicinity of {s^)k will be taken into account in the approximate 
evaluation of the integral, means that the integral over can be replaced by integration over a new variable A, with 
±s° replaced by A in the integrand. This new integral is dominated by the region near Aj. = Ks*^)*:], determined by 



= 



dDkjX) 
dX ' 



(23) 



5 



in which Dk{X) is given by Eq. (|l9| ) with = ^A, L-y/ — . Expanding -Dfc(A) to second order about Aj., the 
rest of the argument of the exponential to first order, the rest of the integrand to zeroth order, performing the A 
integration, and squaring the amphtude yields the result 



'Is+Fs 

n 



•0k. (x5,a:g;) 



4(Vs)^(Vd)^ 

41 19 



1/2 



X exp [-iXkiVD - A) + isk{Xk)\yD -xs| - CkiXk^xs.yo) - DkiX^)] 



X Mn 



Afe - Sfc(Afe)cr • L 



Ml 



(24) 



where 



2 rfA2 



(25) 



Study of the explicit expression for shows that it is essentially the sum of the squares of the time and length scales 
of the wave packet overlaps in the source and detector, i.e. 



il ^ [Tsf + {Lsf + {Td? + {Ld)\ 

(in rather obvious notation). This is particularly transparent in the limit of relativistic neutrinos. 
The factor exp[— Cfe(Aj,, xs, j/d)], with 



(26) 



1 



Ck{Kk,xs,yD) = ^ 



D -^Sl 



1 

Vk 



Wd - xsl 



(27) 



suppresses contributions from neutrinos that do not follow a classical spacetime trajectory between the production 
event at (a;g,X5) and the detection event at (jj%,yD) [the neutrino velocity is given by Vk = Sfc(Aj,)/Aj,].^ 



IV. MACROSCOPIC EVENT RATE 

To make contact with experiments it is necessary to magnify the probability of Eq. ( p^ up to macroscopic 
scales. For this purpose, the normalization in Eq. (^ suggests that (for example) |^/;k. (x5,Xg)| be interpreted as 
the (per particle) volume density of particles with momentum kj at position xs and time a;^, where the last two 
quantities are now thought of as macroscopic spacetime variables. Employing the usual statistical methods for free 
particles, these particle densities are taken to be [dkj/(27r)^] /(kj, xg, a;^) for initial state particles (where / is the 
phase space density) and [dkj/(27r)'^] for final state particles. In connection with the (now macroscopic) variables 
Xs, x'g, Yd, and y^, one factor of (VsVd) is interpreted as dxs dx'g dyn dy'jj. At the macroscopic level, a sum over 
external particle spins is performed; the average over initial spins is accounted for by leaving the spin degeneracy 
out of the phase space distribution functions /. The expected number of events detected from neutrino interactions 
with the [dp/(27r)'^]/(p, y£), y^) particles of momentum p in detector volume dyo during time dy'jj resulting in 
final state detector particles of momentum {Pj,}i arising from neutrinos produced from the interaction of the set of 

{[c?kj/(27r)'^]/(kj, X5, Xg)} detector particles with momenta {kj} in source volume dxs during time dx^ resulting in 
final state source particles of momentum {k^/}, is 



^The above procedure in which a "neutrino energy" At. is determined only from the minimum of Dfc(A) implicitly assumes 
that the phase ^\{y% — x%) + Sk{\)\yD — x^] in Eq. is essentially stationary over the range of A for which exp[— _Dfc(A)] is 
appreciably nonzero. It is easy to see that the quantity exp[— Cfc(A^, xg, yo)] enforces this very condition, so that the procedure 
is self-consistent. 
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dN = dxs dx% djD dy% dK{xs) dK' dPiyo) dP' 



71-'' I Yd 



X5l 



■E 

spins 



1/2 



X Mn 



where the notation 



dKixs) 



dK' 



Ml 



Is 

n 



dk, 



Fs 

n 



;27r)3 (2i?k^ 

,A (27r)3 (2i?kJ 
dp 



(27r)3 2Ep 



f (p,yi5,y?)) , 



dp ., 



,v (27r)3(2i?p^,) 



(28) 

(29) 

(30) 
(31) 

(32) 



has been introduced for the phase space factors. 

While virtually all neutrino experiments record data that sums over contributions from all initial momenta in the 
source and detector, all source final momenta, and all source emission times, some experiments — such as those with 
water Cerenkov detectors — record the detector event time and (at least some) detector final state particle momenta. 
Hence it would not be correct to integrate over these last quantities. Also, dividing by dy'^ gives an expected event 
rate as a function of detector time y^,. 

In the integration over all source times Xg, the formalism automatically "knows" that neutrinos emitted at macro- 
scopically different source times are not allowed to interfere coherently. Each term in the squared sum of the form 
J2k — J2k h{k)h*{k') in Eq. (|2^ ) has a factor exp[— 7^^'], where 

Tkk' = -i{kk -Xk'){VD - A) -Ckikk^xs^yo) - Ck'{\y .xg^yo)- (33) 

Terms with k ^ k' represent quantum interference terms. For a given time of detection y^, exp[— Cfc(Aj., xs, yi))] 
and exp[— Cfc' (Aj,/ , xs, J/d)] tend to pick out different emission times for k ^ k'. If the difference in emission times is 
greater than the width of exp[— C], the interference is suppressed. 

The gradual loss of coherence can be expressed quantitatively. The leading contribution to terms with k ^ k' comes 
from the interval near the average emission time 

(4)fcfc' lyu-xsl, (34) 

Hkk' 

where Ck{X^.,xs,yD) + C'fe' (Aj,/ , xg, yc) has a minimum. The "average velocity" Wj.^., is given by 



V-kk' 



[elvk + il,vk') ■ 

The portion of the argument of the exponential that depends on — that is, Tkk' of Eq. ( ^3| ) — can be expressed 



(35) 



^ _ -^(Afc-AfcO , 

■ikk' — IYD 

V-kk' 



xsl 



lYD - X5| 



i^k-h'mel {^1 + 4') 



X% — {Xg)kk' — 



2iel£l,iXk-Xk') 



(36) 



The second term yields an exponential fall off with \yo — x, 



(and hence different velocities). Setting Eqs. (|3^ ) and (|6 

external momenta, and dividing by dy'^ gives the expected event rate in the detector at time y^. 



in interference between neutrinos with different masses 
into Eq. (p8|), integrating over x^ and the unobserved 



^Unlike the case of Eq. ( |l8| ) in which there was a self-consistent way to assume a stationary phase, here it is necessary to 
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V. RELATIVISTIC LIMIT 



Here the event rate in the hmit of relativistic neutrinos wiU be exhibited. The zeroth order neutrino energy A and 
coherence width £ are given by Eqs. ( psf ) and ( p5| ) respectively, with mj. = 0. To first order, 

Afe = A + 5Afe, 

2 

ID 

Sfc(Afc) = A + <5A,-^, 



Vk 



4 



Ikk' 



1 _ Illfc 

2A2' 

e + s£i 

_ {ml + ml,) 
4A2 ' 



(37) 



where the expUcit forms for SX^, and 5£l are determined by Eqs. (^3|) and ( |25| ) respectively. It will be assumed that 
m| can be neglected everywhere except when appearing with the macroscopic distance lyp — X5I in the argument of 
the exponential, which magnifies its impact. This means that the third term in Eq. ( |36| ) can be neglected, and that 

DkiXk) « D{X), (38) 



where -D(A) is given by Eq. (|19| ) with S^k replaced by ^ = I A, AL) . It also means that 



Afe - Sfe(A;,)cr • L 



Ml «M2 



A - Act ■ L 



Ml 



= [-M5({k},A)] MDi{p},X 



(39) 



where Ms ^^nd Md are the S function- free matrix elements that would appear in the plane wave S'-matrices [see Eq. 
(|l|)] describing the source and detector reactions with a massless neutrino of momentum AL, with the other particles 
having momenta {k} (source) and {p} (detector) |8j. The event rate at detector time y^) is 



SVsVl 



E l-^s({k},A)| 



spins spins 

x^UakUpkUlk-Unk' exp 

k,k' 



\2P) 



{ml-ml,)\yD -^s\ {ml - ml,Y\yD - 



121 



2A 



32A*^2 



(40) 



It can be shown that 



-2D{\) 



(41) 



and consistency with the earlier approximate evaluation of the (or A) integration means that |Als ({k}, A)|^, 
A1_D ( {p}, A ) , and the factors summed over fc, k' can be taken inside this integral as functions of iJq rather than A. 



include the phase in completing the square for Xg. It is this which gives rise to the second term in Eq. (^), which causes a 
loss of coherence with increasing ly^ — xgp. It can be shown that the exponential fall off in \yo — x^p also ensures that the 
phase — A(i/^ — x^g) + Sk(X)\yD — xs| in Eq. ( |l^ ) remains stationary over the range of A for which exp[— _Dfc(A)] is appreciably 
nonzero, for x% determined by the average velocity v_f.y [see Eq. (|34|)]. It will be recalled that the stationarity of this phase 
ensures that the minimum of -Dfc(A) dominates the integral (i.e., the "neutrino energy" becomes a meaningful concept). 
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In addition, if the phase space densities change little with energy variations and momentum variations of order i ^ ^ 
then the leading contribution to Eq. (^ ) is the same as if the replacement 



Dct [(M^s^)^.]y'Dct {{W^^)^. 
■ 5^ i-ks + q) 5^ {PD - q) 



: i-ks + q) 5^ {pD - q) 



VsVl 



(42) 



had been made, where q — (£^q,i?qL). Hence the leading contribution to the macroscopic event rate in the detector 
at time y^, can be expressed 



dp 



n 



(27r)3 



yV,-|yi3-xs| 



(43) 



where the single particle event rate is 



dr({k},{p},xs,yc) = J dE. 



dr({k},i?q) 



Yd - xspdfiq dEc^ 
X [Pmix (£^q, xs, Yd)] da (^{p}, i^q^ 



(44) 



In Eq. (El 



dEa 



rfr({k},i?q) 



Iyd - xsp dflqdEq 



i?q di?q 



Fs 

|yz,-xs|2(27r)3(2i?q) 1^ ^9R. ^ 11 
X ^ |Xs({k},i?q)|' (2^)V(-fc5 + 9) 

spins 



Is ^ 



dk,/ 



(27r)3 (2i?k,,; 



(45) 



is the flux of neutrinos of energy Eq at position yo due to an interaction at X5 , as would be computed with standard 
plane wave methods; 



da 



({P},£^q) 



1 



(2Sq) {2Ep 



n 



dp ., 
—J 



(2^)3 (2£;pJ 



Y,\MD({V}.E^ {2n)H'{pD-q) 



(46) 



is the cross section for a massless neutrino interaction in the detector (assuming nonrelativistic initial state detector 
particle momentum p, so that the M0ller velocity is equal to 1), and 



Pniix (i;q,Xs,yz5) = ^ Ua^kUpkKk'UfSk' 
k,k' 

.{ml-ml,)\yD-^s\ {ml~mlf\yD-^s? 



X exp 



2En 



32E^£^ 



(47) 



is the flavor mixing (or "oscillation" ) probability. 

Except for two differences, Equations (^3|)-(^7|) are just what one would write down for a macroscopic event rate 
using the naive QM model of the neutrino flavor mixing process, together with elementary considerations for the 
production flux and detection cross sections. The first difference is one that could also have been put in by hand. 
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namely, the causal time delay between emission and detection.]] The second difference is the damping of coherence at 
very large distances, discussed earlier in this section. 



The calculations presented here concern the description of neutrino flavor mixing as a quantum field theoretic 
process, with the neutrinos as virtual particles connecting the "on shell" external particles involved in the neutrino 
production and detection reactions. This framework provides a more realistic description of and deeper physical insight 
into the flavor mixing process than the naive quantum mechanical model. Development beyond previous works has 
been sought in this study by considering general production and detection processes and wave packet functional forms, 
avoiding the unrealistic assumption of microscopic stationarity, leaving the time of the detection event (but not the 
emission event) an observable quantity, and making a complete connection to fully normalized event rates. The final 
results are given in Eqs. (^)-(^). Note that Eq. ( ^7| ) has a damping factor for interference terms in addition to 
the usual oscillatory factor; this should formally be considered part of the "oscillation probability," which (without 
ambiguity, in the present fully normalized treatment) includes everything outside of the neutrino emission flux and 
detection cross section. 

The free external particle wave packet picture is convenient for a number of reasons. Unlike descriptions involving 
bound states in the source, it is a suitable description for astrophysical neutrino sources such as the Sun (modulo the 
Coulomb repulsion of reacting nuclei, which can be suitably included) or supernovae. The lack of stationarity arises 
naturally due to the finite duration of the wave packets' overlap. The matter of normalization is simple since free 
particle states are employed. Furthermore, the dynamics are already built into the ^-matrix, making the description 
of neutrino oscillations a matter of working out the kinematics. 

This S-matrix framework could be generalized without much difhculty to include bound states for some of the 
source and detector particles. (The remaining particles would still be considered free particle wave packets. In this 
framework, the coherence times of the source and detector — which ultimately result from complicated microscopic 
many-body physics not considered here — can be considered as input parameters which ultimately manifest themselves 
in the finite free particle wave packet sizes.) The analogue of the S'-matrix would be the amplitude for particular 
plane wave states to interact with particular bound states. A superposition of such amplitudes over several plane 
wave momenta — in order to create time-dependent wave packets for the external free particles — would constitute 
the amplitude for the overall neutrino production/propagation/detection process, with the neutrino production and 
detection localized in space and time. In going to the macroscopic rates, the square of bound state wave functions 
would be replaced by a sum over the phase space distribution of the relevant bound state quantum numbers. 

While not new to this study, three basic insights into the neutrino flavor mixing process are listed here for com- 
pleteness. First, (1) an "oscillation probability" independent of the details of production and detection can only be 
defined in the relativstic limit. This limit allows the neutrinos to become effectively on shell (i.e. massless) as far 
as production and detection are concerned. Assuming chiral interactions, the relativistic limit also causes only one 
neutrino spin to contribute, so that the overall squared amplitude |A^p can factorize into separate production and 
detection squared amplitudes |A^sP and |A^Dp. This process has been shown in detail in this paper, culminating in 



The second insight is a condition on (2) the maximum size of the external particle coordinate space wave function 
overlap in the source (Ls) and detector (Ld) that allows neutrino states of different mass to interfere coherently: 



When expressed in terms of the "oscillation length" this condition is intuitively obvious. Its necessity can be seen 
mathematically in Eqs. (|l^) and (p4|). From these equations one can see that the "neutrino energy" is determined 
by a compromise between the degrees to which energy and momentum are conserved in the source and detector. 
However, the tendency towards energy conservation generally has a greater impact, i.e. the timelike components 
(W^^)oo {TsY and (VF^^)oo {TdY (where T5, Tjj are the time scales of the wave packet overlaps in the source 



''In setting x% — i/^ — \y d — y^sl in Eq. ([l3|), it has been assumed that the phase space densities / vary on time scales slower 
than 8 x 10~^^ s [(m| + my) I eV^](MeV / -Eq)(jy_D — xs] / km), for all k and A:', in order that the phase space densities 
could be taken out of the sum. 



VI. DISCUSSION 



Eqs. (g|)-(g3). 




(48) 
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and detector) are larger than the spacehke {Wg^)ii ~ [LsY and {W^^)ii ^ {Ld)'^ since the external particles travel 
slower than the speed of light. To the extent that the energy is more well determined, there must be a greater spread 
in momentum in order for interference to occur, which is why the condition expressed above is couched in terms of 
the spatial (as opposed to temporal) spread of the wave packet overlap. 

The finite duration of the production and detection processes leads to the third insight. Contributions to the 
amplitude from neutrinos that deviate from a classical spacetime trajectory are exponentially suppressed [see Eqs. 
( p^ ) and (p^)]. This leads to (3) an upper limit on the number of observable oscillations in space: 

N^^^Jl2_^<,-^E.i^^ (49) 



[see Eq. (|7|)]. Here the "coherence width" i ~ ^J{TsY + {Tof + {Lsf + {LdY [sec Eqs. (|2|) and (|6|)], and 
the detector resolution AiJ,y = Ap^j has been taken as a crude estimate of this quantity. Thus many oscillations in 
space should be visible before decoherence sets in as the spatial and temporal^ resolution of the detector begins to 
distinguish the separating neutrino mass eigenstates. 

It is true that the three insights above can, to some extent, be achieved without the elaborate machinery presented 
here. As far as insight (1) goes, common knowledge that weak interactions are V — A makes the irrelevance of the 
neutrino spin degree of freedom in the relativistic limit somewhat obvious. The relativistic limit also makes the 
notion of "real" flavor eigenstates reasonable (zero mass is on-shell). One can then adopt the simplified quantum 
mechanical picture. In connection with this simplified quantum mechanical picture, condition (2) can be argued from 
the uncertainty principle, AxAp ^ 1, where Ax corresponds to the oscillation length Lose and Ap corresponds to the 
inverse source/detector sizes [Q. Condition (3) follows from noting that a real source (detector) will have some 

finite linewidth (resolution). In order that interference terms not wash out when binned over this energy range Ai?, it is 
necessary (e.g., [|l^]) that the variation in the oscillation phase, A[(m^ — TO^,)|y£) — xs|/i?] ~ ([yc — ys|/iosc)(Ai?/i?), 
be smaller than 27r, which is essentially condition (3). 

Even if these arguments can be made in some fashion in connection with the simplistic formalism, the whole picture 
lends itself to conceptual difficulties at some level Furthermore, the conditions (l)-(3) must be invoked from prin- 
ciples outside the formalism itself. In contrast, the QFT description of the neutrino production/propagation/detection 
presented here exhibits all of these conditions in a natural, self-contained, and physically unambiguous manner, and 
describes the transition to the failure of these conditions in a quantitative way. 
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